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I Abstract 

Fluctuations of the interface between coexisting colloidal fluid phases have 
O been measured with confocal microscopy. Due to a very low surface tension, 

I— —I the thermal motions of the interface are so slow, that a record can be made of 

^ the positions of the interface. The theory of the interfacial height fluctuations 

^ is developed. For a host of correlation functions, the experimental data are 

compared with the theoretical expressions. The agreement between theory and 
experiment is remarkably good. 

O 1 Introduction 

00 

^ The study of interfaces has a long and interesting history. In 1894 van der Waals [1] 

. ^ proposed an interface theory, which leads to a fiat interface with a density profile in the 

^ direction of gravity. This result is sometimes referred to as the intrinsic interface. Al- 

ready von Schmoluchowski [2] realized that the thermal motion of the molecules induces 
height fluctuations in the interface. These motions have been called capillary waves, 
since they derive from an interplay of gravity and surface tension, like capillary rise. 
The fluctuations were first treated theoretically and experimentally by Mandelstam |3]. 
He pointed out that the interface width diverges due to the short wavelength capillary 
waves. This fact was rediscovered by Buff, Lovett and Stillinger |1] fifty years later, 
after which it obtained a prominent place in the discussion of interfaces. Weeks [5] 
later pointed out that the notion of capillary waves only applies to wave lengths larger 
than the fluid correlation length, which is of the order of the interparticle distance. 

The experimental study of interfaces was undertaken by Raman [6| and Vrij [7] 
with light scattering and starting with Braslau et al. [H] by X-ray scattering. Although 
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Figure 1: Schematic view of Confocal Microscopy. The confocal microscope thin focal 
planes of approximately 500 nm thickness can be imaged. This enables the investigation 
of local phenomena such as height fluctuations. 



scattering on interfaces is most valuable, it always yields global information on the 
fluctuations, while a photographic inspection gives local information. However, the 
wave lengths and heights involved in the capillary waves of molecular fluids are way out 
of the reach of detection by photographic methods. The visual inspection of capillary 
waves initially remained restricted to computer simulations of interfaces in molecular 
systems [9j. 

The field obtained another dimension by recent experiments of Aarts et al. p[0] in 
which they obtained pictures of fluctuating colloidal interfaces. The key is that, by 
lowering the surface tension to the nN/m range, the characteristic length and time 
scale of the fluctuations become accessible by confocal microscopy. This opened up 
the possibility to follow in detail the motion of the height of the interface and to do 
a statistical analysis of its temporal and spatial behavior. Of course the method has 
its inherent restrictions. Just as in ordinary movie recording, the pixels have a finite 
distance and the snapshots have to be taken at finite time intervals. For colloidal 
interfaces this interval can be made much smaller than the intrinsic time scale of the 
motions. Thus the Brownian character of the motion could be demonstrated ad oculos. 

In the confocal microscopy a two-dimensional section is inspected perpendicular to 
the interface and the density profile between the two phases is observed. A schematic 
picture of the experiment is shown in Fig. [Tj A very precise location of the interface 
can be obtained by fitting the intensity with a van der Waals-like profile: I{z, x) = 
a + 6tanh([z — h{x)]/c), where z is the direction perpendicular to the interface and x 
a coordinate along the interface. In the upper phase the density approaches a value 
corresponding to a + 6 and in the lower phase to a — b, while c measures the intrinsic 
width of the interface. Thus at every snapshot a function h{x) follows and the sequence 
of snapshots gives the function h{x,t). This is a practical separation of the particle 
motions, which lead at short scales to the intrinsic interface and the particle motions 
which drive the long wavelength capillary waves. One might think that, if the time 
interval of the snapshots is sufficiently small with respect to the characteristic time 
scale of variation in h{x,t), one can analyze h{x,t) as a continuous function of the 
time, like a movie gives the impression of a continuous motions, while it is a succession 
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of snapshots. In a previous short report [l_3j on these experiments we have pointed 
out that the statistics remains dependent on the time interval, due to the Brownian 
character of the motion. 

We have performed confocal microscopy measurements on phase separated colloid- 
polymer mixtures. The colloids are 69 nm radius fiuorescently-labeled polymethyl- 
metacrylate particles, suspended in cis/trans decalin, with polystyrene (estimated ra- 
dius of gyration = 42 nm [TT]) added as depletant polymer. Due to a depletion induced 
attraction these mixtures phase separate at sufficiently high colloid and polymer volume 
fractions and a proper colloid to polymer aspect ratio into a colloid-rich/polymer-poor 
(colloidal liquid) and a colloid-poor /polymer-rich (colloidal gas) phase [12]. Here the 
polymer concentration acts as an inverse temperature. By diluting several phase sep- 
arating samples with its solvent decalin, the phase diagram was constructed. With a 
Nikon E400 microscope equipped with a Nikon CI confocal scanhead, series of 10 000 
snapshots of the interface were recorded at constant intervals At of 0.45 s and 0.50 s 
of two statepoints to be denoted as II and IV. The pixels are separated by a distance 
Ax = 156nm and a single scan takes approximately 0.25 s to complete. 

The set-up of the paper is as follows. We start out by discussing the spatial be- 
havior of the data of a single time frame, which requires only equilibrium statistics. 
The correlation functions and the statistics of hills and valleys in the interface are 
determined and compared to the theory. 

Then we identify the set of interface modes via the fourier decomposition 

/i(x,t) = ^/ik(t)exp(ik-x). (1) 

k 

The modes are overdamped in the relevant regime and follow from the macroscopic 
interface-dynamics. The motion obeys not only the macroscopic equations but is also 
influenced by noise. We introduce thermal noise through the Langevin equation and 
calculate the essential height-height correlation function {h{0, 0)h{'x, t)). Via the equiv- 
alent Fokker-Planck equation the probabilities on sequences ( "histories" ) of snapshots 
are determined. The analysis of the distributions of "hills" and "valleys" in the time 
domain with respect to a level h is similar to the spatial behavior. A special concern is 
the dependence of the residence time and the waiting time on the used time-interval. 
The paper closes with a discussion of the main results. 

2 Equal time Correlations 

The function /i(x, t) provides a mathematical division between the two coexisting 
phases which form the interface. The interface is of the solid-on-solid type since so- 
called overhangs, well-known in lattice theory, are excluded by construction, as to every 
value of the horizontal coordinate x and time t one unique height /i(x, t) is associated. 
The basic function is the height-height correlation function. Due to translational in- 
variance the modes k are independent and thus the correlation function in space has 
the fourier decomposition 

{h{0, 0)/i(x, 0)) = ^^P(^k ■ x) (2) 

k 

The brackets denote equilibrium averages and /ik is the amplitude of the k-th mode. 
The distribution of the /ik follows from the Boltzmann factor involving the energy of a 
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Figure 2: The height distribution for statepoint II and IV as found experimentally. 
The lines are gaussian fits to the data. 



deformation of the interface, which is given by the drumhead model 

= \j [^Pgh\-^) + 7(V/i(x))2] _ 

Here Ap is the density difference between the coexisting phases, 7 is the surface tension 
and g is the gravitational acceleration. The first term gives the gravitational potential 
energy and the second term the increase of the interfacial energy due to curvature. 
Expressed in terms of the amplitudes /ik it reads 

mh}) = ^Y.^^P3 + ^ (4) 

k 

where is the area of the interface. Since Q is quadratic in the amphtudes /ik, it 
implies a gaussian distribution for the /ik 



with the average 

With the distribution (|5| of the h^, we can calculate the distribution of the heights 
h, which becomes also a gaussian 



4 



with the mean square height (fi^) given by 



^ is the capillary length defined as 



The integral has been given an upper bound k^ax — iT^jd with d the diameter of the 
particles and a lower bound kmin — 2'k/L due to the finite size of the interface. The 
lower bound can be set equal to for all practical purposes, but the upper bound 
is essential for the convergence of the integral. Cutting off the capillary waves at the 
short-wavelength side is the poor man's way to handle the otherwise diverging interface 
width (/i^). There are two options to determine {h'^)- The first follows from a fit to 
Pgq{h), which is shown in Fig. |2| The second is a direct evaluation of {h'^) from 
the recorded data. The latter always gives a 1-3 % larger value, which we attribute 
to optical artifacts due to confocal slicing. So we are inclined to prefer the former 
value which amounts to {h"^) = 0.219 for statepoint II and 0.336 {fim)'^ for statepoint 
IV. Then equation ([s]) can be used to estimate the upper cut-off. On the basis of a 
determination of 7 (see below) one finds values around k = k^ax^ — 45, but this value 
is rather sensitive to small variations in 7: a variation of 7 by 10-15% results in a shift 
in K by a factor 2. 

The correlation function (/i(0, 0)/i(x, 0)) is of course also measurable. In the ap- 
pendix we discuss the integral (|2]); here we give the result with the cut-off sent to 
00 

(Mo,o)Mx,o)) = ^iro(x/0. (10) 

The divergence for a; ^ of the modified Bessel function Kq corresponds to the 
divergence of the interface width without a cut-off. A fit of the correlation function to 
the Besselfunction (with a slight modification due to the cut-off) is shown in Fig. |3] 
Apart from a few initial points the function fits quite well. We find fitting parameters 
^ = 8.0/im for statepoint II and ^ = 6.1/im for statepoint IV. The values for 7 turn 
out to be 58nN/m viz 21nN/m for statepoint II viz. IV. 

Fourier transforming the correlation function back to the wavenumber domain 
should lead to the expression (|6| as function of k. However, an inverse fourier transform 
requires accurate data for a large domain and the correlation function is unreliable for 
large distances (not shown in Fig. [3]). This prevents a direct check of the drumhead 
hamiltonian. 



3 Multiple Correlation Functions 

As the data are stored for all sampled positions we can determine the probability 
density 



G'n(/ii, Xi; ■ ■ ■ ; hn, x„) — 

(<5(/l(Xi, 0) - /li) ■ ■ • 5(/l(x„, 0) - hn)), (11) 

which gives the joint probability that the interface at position xi has the height hi and 
subsequently at position Xj the height hi etc.. A straightforward evaluation of (11) 
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Figure 3: The spatial correlation function (/;,(a;, 0)/i(0, 0) for the statepoints II and IV 



fitted to the expression ( 10 ) using a cut-off. The lines are fits to the experimental data. 



proceeds via writing the 5 functions as a fourier integral and then expressing /;,(x, 0) in 
terms of the amplitudes /ik- As all integrals are over a quadratic form in the exponent 
the result is a gaussian in the hi. So one gets an the expression of the type 



G'„(/ii, Xi; 



det g 



-1 



(27r(/i2))' 



1/2 



exp 



1 hihj 

2 Z^^M Jf^ 



In this notation the matrix Qij is the correlation matrix 

9i,j = gil^i - Xj|,0), 
which turns out to be the equal-time height-height correlation function 

^(|x-x'|,t-t') = (Mx,t)Mx',t'))/(/^')- 



A shortcut to the result (12) is to evaluate the following integral in two ways 



dhi--- dhn hihj G{hi, Xi; ■ ■ ■ ; /i„, x„) 



(12) 



(13) 



(14) 



{h')g,j. (15) 



The first uses the definition (11) and obviously leads to the right hand side of (15). 



The second way uses the expression (12). Then one has to diagonalize the quadratic 



form in the exponent and the integration over the eigendirections also leads to the right 



hand side of (12), which shows that (13) is correct 
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In (14) and accordingly in the result (12), we have factored out (h?) because we 



want to use it as a scale for the heights. (12) shows that the correlation function g 



dictates the behavior of the multiple correlation functions. 

To give an impression on the behavior of the G"s we first consider a few small 
values of n. For a single position {n = 1) the value g{0, 0) = 1 and (12) reduces to the 
equilibrium height distribution ([s]) 



(16) 



which is shown in Fig. [2] From this probability we derive the important probabilities 
q~^{h) to find a height above h and q~{h) for finding a height below h. They are given 
by the expressions 



g+(/i) = / Pe,ihi)dh, 
'h 



q-{h) 



Pe,{h)dh. (17) 



In integrals like (17), one changes of course to the combination hi/ (/i^)^/^ as integration 
variable, such that the q'^{h) become functions of the scaled variable h/ih?)^^"^. The 



result of the integration in (17) is an error function in this parameter. From now on 



we work with these reduced heights. 
The probability density G for n - 



2 is still sufficiently simple to make it explicit. 



The matrix gi^j and its inverse g^ j read 



9i,j 



1 9l,2 
92,1 1 



(1^ 



and 



9^ 



1 - 9I2 



1 

'92,1 



So G2 follows from the general definition as 

1 



G'2(/ii,0;/i2,x) 



27r[l - C/2]l/2 



exp 



'9l,2 

1 



hi - 2ghih2 + hj 



(19) 



(20) 



with g a shorthand for gi^2 = g{x, 0) and x the distance of sampling. Note that this 
expression is symmetric in the entries hi and /12 and that dependence only enters 
through g = g{x, 0). 



4 Statistics of Sequences 

The probability densities ([TT| are measurable, but the statistics becomes poor when 
too much entries are taken. Therefore integrated probabilities are more accessible. 
For what follows it is interesting to study the probability that a sequence in space of 
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Figure 4: The function q^^{h, g) as function of x through g = g{x, 0), for some values 
of h (in units (/i^)^/^). The drawn hues are calculated and the points are the measured 
values. 



precisely n successive values occurs of the heights above the level h. According to the 
theory it is given by the ratio of two integrals 

Pi{-h)=q~^''^^'ih)/q^-{h). (21) 

In this notation the superscript prescribes the integration domain. The numerator of 



(21) reads 

/h POO 
dho / dhi-- 
-oo J h 



/■oo rh 

/ dhn dhn+iGn+2{ho,- ■ ■ ,hn+i)- (22) 

J h J — oo 



The integral over the first variable ho guarantees that the sequence starts below level 
h, the next n integrations select points above the level h and the sequence ends with 



hn+i below level h. So the numerator in (21) selects the hills of precisely n consecutive 



values of the height above level h. We have omitted in Gn+2 the position arguments 
since it is understood that points are equidistant. So the sequence of values g{mAx, 0) 
enters, with < m < n + 2. The denominator is the integral 



g+-(/i) = / dh, dh2G2{h,h2). (23) 

J h J — oo 

It counts the number of hills since each hill is followed by a transition from above to 



below the level h. The denominator serves as a normalizing factor. Summing (21 ) over 



n (from 1 to oo) gives the total number of hills above h and as this equals the number 



of crossings, we see that the distribution (21) is normalized. 
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Figure 5: The spatial Pn{h) as function of h for a number of n for statepoint IV. The 
drawn hues are the theoretical values. 



The average length x^W of a sequence is defined as 



(24) 



n=l 



We also look to sequences below the height h. They are given by the probability 



Pn{h) which follows from a similar definition as (21), with + and — interchanged. The 
up-down symmetry of the problem yields the relation 



Pnih) =p^{-h). 
The average length x {h) of a stretch below h similarly equals 

n=l 



(25) 



(26) 



Inserting the expression (21) into the definition (24) for x^{h), the numerator in (22) 



is multiplied by the number of values larger than h. Summation over all n leads to the 



average number of points above the level h, which is given by the integral (17). Thus 
we arrive at the relations 

X^ih) = q^{h)/q^~ih). (27) 

The remarkable point about these relations is that, although the probabilities p^{h) 
are given by multiple integrals, the averages x^(^) result from simple integrals. The 
q^{h) are errorfunctions and is a two- fold integral involving the function G2, thus 
containing only the value 5'(Ax, 0). 



A trivial result from (27) is that the ratio X~^{h)/x~{h) is the same as the ratio 
q^{h)/q~{h). Both give the ratio of the total number of points above and below the 
level h. As the q^{h) and q~{h) add up to 1, a more intriguing result follows for the 
sum 

X^ih) + x-ih)= ^ \ (28) 

{h,g) 
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Deliberately we have given q^^ also the argument g which incorporates the spatial 
dependence on x. Tacitly we have assumed that this distance is the sampling distance 
Ax. But nothing prevents us from taking a multiple n of As. Then g will refer to 
5f(nAa;,0). In Fig. |4]we have plotted the experimental values of q^^{h,g) for various 
values of g{nAx,0), which we take as parameter on the horizontal axis. The curves 
are the calculated values of q~^~{h,g). We have not found a closed expression for q~^~ 
in terms of known functions, but a number of limits are explicitly obtainable. The g 
dependence is exemplified by the case h = 0, which reads 

(0,9) = 7;-- arctan . (29) 

2 vr \i--9j 

The h dependence is by and large controlled by the limiting behavior 

q^-{h,0) = q+{h)q-{h), 

./T^ r 1 

(30) 



q [h, g-^1}^ T=- exp 

TTV 2 



Apart from the averages also the individual (h) can be measured and compared 
with the theoretical expressions (21). In Fig. [sjwe show p^(/i) for statepoint IV as a 



function of h for a number of n. The theory requires the evaluation of multiple integrals 



(22) which can be carried out by Monte-Carlo integration. The best procedure is to 
generate a distribution according to the gaussian integrand and then reject the points 
that fall outside the integration domain. This technique also applies to correlation 
functions for other histories with another integration domain. 

The comparison with the theory has been carried out for n up to 8. For higher 
values of n the integration becomes a bit lengthy for good statistics. The agreement 
between theory and experiment is good, but typically there are deviations for larger 
negative values of h, where the experimental points are systematically lower than the 
theorical prediction. 



5 The Dynamic Interface Modes 

As the preceding sections show, the interface fluctuations have a rich spatial structure. 
So it is an interesting question how this compares with the interface fluctuations in 
time. The temporal development of the interface is determined by the macroscopic 
equations for the interface modes as well as by the influence of thermal noise. In this 
section we briefly discuss the interface modes and in the next section we treat the noise. 

The problem of the interface modes modes has been addressed by Jeng et al. [Hj, 
who have made an extensive study of the interface modes in the various regimes dis- 
tinguished by the relative strength of viscosity and surface tension. The modes are 
overdamped for our experimental conditions and decay as 

hk{t) = hkexp{-ujkt), (31) 

with a rate 

cok = ^m)-'+m, (32) 



where the capillary time tc is given by 



(33) 
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Here rj and 1]' are the viscosities of the lower and upper fluid. A few remarks on (31) 
are worth making: 



The dispersion relation Uk is in general rather complicated. The simplification 
(32) derives from the approximation puk <C rik"^ which is very well fulfilled for 



colloidal interfaces with extremely low surface tension. The approximation is 
controlled by the number 



\ 9 J iv + V' 



(34) 



which is the ratio of the capillary length ^ to the viscous length = [rj + 
ri'y /'-fAp. It is very small, 10~^, for colloidal interfaces, while it is very large for 
e.g. water (10^). 

The spectrum has a slowest mode with wavelength ^ and decay rate tc, in con- 
trast to the capillary waves of molecular fluids, where the modes become slower 
the longer the wavelength. This mode starts to dominate the behavior of the 
correlations for long times. 



6 The Langevin Equation 

The thermal influences can be incorporated by a fluctuating force F]^{t) on mode k in 
the Langevin equation [T5] 

-^ = -oOkh^ + F^{t). (35) 
ot 

The first term on the right hand side is the systematic damping force, which by itself 
would lead to an exponential decay of mode k. The random force F]^{t) has zero average 
and is assumed to be ^-correlated in time (white noise) 

(Fk(t)Fk'(t')) = W,ork5(t-t'), (36) 
where can be found from the fluctuation-dissipation theorem 

M). (37) 



_ |2\ 



The Langevin equation (35) assumes that the slow capillary waves form a complete 



set to characterize the motion of the interface. is linked in (37) to the equilibrium 
average of the amplitudes /ik, which is given by 

With the Langevin equation all time-dependent correlation functions can be calcu- 
lated. In particular the height-height correlation function follows as 

{h{O,0)h{^,t)) = {h')g{x,t) = 

^(l/ikH exp(ik ■ X - t^fct), (38) 



with (|/ikP) given by ^ and (/i^) by (jSj). Note that for this correlation function the 
influence of the fluctuating force -Fk(t) averages out such that it is depends only on the 
macroscopic dynamics of the interface. It involves, apart from the decay rate Uk only 
the thermal average (|/ikP)- Some properties of the integral yielding this function are 
listed in the Appendix. 
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Figure 6: The correlation function (/i(0, t)/i(0, 0)). The points are experimental values 



and the lines curves according to (40) using a cut-off. 



The first point is the determination of 7 and tc from the data. We represent 
{h{O,0)h{O,t)) as 

(/i(0, 0)h{0, t)) = ^Hit/t,, «:), (39) 



27r7 



Here again k = k^axi- In the Appendix we prove that 
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(40) 



with Kq the modified Besselfunction of order 0. For t > tc/n, the function H{t/tc, k) 
is well represented by Ko{t/tc). Since k is of the order 40 to 50 (see Section |2]), 



expression (40) suffices for most of the measured points, except of course for the first 



few points near t = 0, where the right hand side of (40) diverges. Leaving them out 
for the moment, we find from a fit for statepoint II: tc = 20s and 7 = Q6nN/m and 
for statepoint IV: tc = 33s and 7 = 22nN/m. Effectively tc acts as a horizontal scale 
parameter and 7 as a vertical shift, tc is mainly determined by the asymptotic behavior, 
while 7 is more sensitive to the initial behavior. 

By adjusting the upper cut-off, the calculated g{0, t) assumes the value 1 for t = 0. 
In Fig. [6] we plot the experimental values of {h^)g{fd,i) together with the theoretically 
calculated curves. 



Finally we mention the initial behavior of g{0,t). From the expansion (62) we 
deduce 

Here one sees that a finite k is essential for this initial behavior. 



7 Probabilities on Histories 



The noise term comes into the picture when we calculate the distribution of the /ik(^)- 
It follows from the Fokker-Planck equation, which is equivalent with the Langevin 
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equation and reads |i5j 



dt 



(42) 



It gives the evolution of the probabihty distribution P{h]f_,t) starting from an initial 



distribution P(/ik,0). The solution [15] of (42) provides the conditional probability of 
the mode h^it), starting with the value /ik(0) 



P(/ik(0)|/ik(t)) 



[27r(|/ikp)(l-e-2-fc*)]i/2 



exp 



|/^k(t)-/^k(0)e-"^T 
2(|/ik|2)(l -e-2-^*) ' 



(43) 



The expression shows that, independent of the value of /ik(0), the distribution asymp- 
totically approaches the equilibrium distribution ([s]). 

For the measurements at different times (and possibly different positions) we need 
the multiple time correlation function 



Gn{hl 1 , ti , ■ ■ ■ , /in, Xn, tj^) 



(5(/l(xi, U) -hi)--- 6{h{^n, h) - hn)), (44) 



giving the probability of a history that the interface is at time ti and position xi at 
a height hi and subsequently at time ti and position Xj at height hi. In order to 
evaluate these correlation functions we have to translate the field /i(x, t) into its fourier 
components /ik- Then we have to use the joint probability on a set of components 
hkiti), - - - h]^{tn), which is given by the product of the equilibrium probability (Isl) for the 



first event at ti and the conditional probabilities (43) for the successive time intervals, 
t2 — ti, ■ ■ ■ ,tn — tn-i- The result of the integration can be derived easily from the 
observation that the h{'Xj,tj) are, as linear combinations of the basic variables /ik, also 
gaussian random variables. So, similar to the derivation of (12), their distribution must 
be of the form 



Gn{hi, - - - ,tr^ 



det g 



-1 \ 1/2 



{27r{h^)y 



exp 



1 i hjhj 

2 ^^''^ (h^) 



(45) 



The matrix gij is the correlation matrix 



9i,j fl'd-'^i 



(46) 



The proof of (45) is exactly the same as that of (12). 



(45) is the main result of the theory for the histories. It relates the probability of a 
history hi, - - - , /i„ on a sequence of snapshots to the height-height correlation function 
g{x,t). The strong point of (45) is that the time-dependent probability densities have 



exactly the same structure as the equal time probabihties, when expressed in the ap- 
propriate gij. Thus the whole analysis given above for the equal time correlations, can 
be taken over for the more general correlations. So we restrict ourselves for the time 
dependent histories to the aspects needing some extra attention. 
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The time dependent probability density G2 reads as (20 ) with g = g(0,t). It can also 
be written as the product of the equilibrium distribution Peq{hi) and the conditional 
probability Gc(^i, 0, 0|/i2, 0, t) that starting at hi one arrives at /i2 at time t later 



a(/ii,0,0|/i2,0,t) ^ 



[27r(l-(72)]l/2 



This expression cannot be seen as the "propagator" for the probability, like ( 43 ) is for 
the fourier components h\^. While the probabilities for the modes k evolve as a Markov 
process, the distribution for h{0,t) does have a memory effect. Only if g{0,t) were 
a pure exponential the spatial process would be Markovian too [I5]. The expression 



(38) shows that 5^(0, t) it is not a pure exponential but a superposition of exponentials. 
For longer times it starts to decay as an exponential when the slowest mode begins to 
dominate, as can be seen in Fig. |6l 



8 Average numbers of hills and valleys 

Consider now a sequence of snapshots, taken with time intervals At. We are again 
interested in the probabilities on the duration of "hills" and "valleys" with respect to a 
level h. To stress the analogy between space and time we use the same notation p^{h) 
for the probabilities to find a stretch of exactly n consecutive values of the interface 
height above/level the level /i, where n now is an index in the time direction. They are 



given by the same integrals as (22) and (23) with gi^j the temporal correlation matrix. 



The mean values are called the residence time T'^{h) (for p^{h)) and the waiting time 
T-{h) {iorp-ih)). 

To check whether the experimental values follow these theoretical predictions, we 



first checked that the ratio T~^{h)/T~{h) is the same q~^{h)/q~{h) in analogy with 27 
It is valid over several orders of magnitude. Only for the very large h deviations occur 
due to poor statistics. In Fig. [7]we now plot again the calculated values of q~^~{h,g) 
as function of the parameter g = g{0,t), for a number of h values. The upper curve 



in Fig. WjfoT = 0, is given by equation (29). In this figure the experimental values 



are plotted as follows. We take as time interval a multiple n of the smallest interval 
At and determine for this sampling rate the r's. This leads to experimental values of 
q~^~{h,g), which we plot in the figure at the value g = g{nAt). The curves for a fixed 
value of h are statepoint independent; the figure shows that this is pretty well the case. 

Finally we plot in Figure |8] the dependence of the r's on h for three choices of 
the time interval. The curves are confusing at first sight, since the values of T^{h) 
are about the same for all three choices. So, if we multiply them with the value of 
the chosen time interval, in order to convert them from numbers to real times, we get 
substantial different times. This indicates that the residence and waiting time depend 
strongly on the measuring process. 

The fact that the smallest chosen time interval leads to the smallest value of the 
residence and waiting time, naturally poses the question what will happen in the limit 
of vanishingly small time interval At. Theoretically it relies on the behavior of the 



correlation function g{t) in the limit t — 0. We presented in (41) the behavior as it 
follows from capillary wave theory. A linear approach of (yf to 1 implies that T~^~{h) 
increases as the inverse power of the square root of At. Then, after multiplying with 
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Figure 7: The function g"*"" as function of t represented hj g = g{0,t), for h = —1 (filled 
squares), h = (circles) and h = 1 (semifilled pentagons). The colours correspond to 
time intervals At = 0.5 seconds (black symbols), 2 At (red symbols) and AAt (green 
symbols). The drawn lines are calculated and the points are the measured values. 



15 




At, 2 At and AAt. 



16 



1t- 




40 aO 120 160 200 

number of intervals n 



Figure 9: The temporal Pn{h) for statepoint IV, for h = —1 (filled squares), h = 
(circles) and h = 1 (semifilled pentagons). The colours correspond to time intervals 
At = 0.5 seconds (black symbols), 2 At (red symbols) and 4At (green symbols) 



At in order to get their values in real time, the residence and waiting time vanish as the 
square root of At. However, the slope of the linear term in (41 ) depends on a molecular 
quantity k, which indicates that wavelengths matter for which the mesoscopic capillary 
waves theory is not designed [5]. One could argue that for molecular times, the cusp 
in g{0,t) is rounded off to a parabola (since it is time reversal invariant). Then this 
parabola would compensate the square root in (30) and the residence and waiting time 
would approach a finite limit. 

Unfortunately this scenario can not be tested experimentally, given the present 
limits on the sampling frequency. However, there is an interesting sampling regime 
beyond our data, for which the capillary wave theory still holds. In Fig. [7] the data 
go up to the value g ~ 0.8. The typical square root decay of the curves for q^~{h,g) 
cannot be tested with our data. A microscope, which is faster by a factor 10, could enter 
this regime where the typical signature of the Brownian character of the fluctuations is 
most significant. They give increasingly larger weight to short living hills and valleys, 
which force the mean values to shrink in a specific way predicted by the presented 
theory. 

The same issue presents itself in the analysis of the data for a single time as function 
of the sampling distance, but in a less severe way. In the appendix it is shown that the 
height-height correlation as function of the distance is a parabola for short distances. 
Thus a finite value of the sequence length and the recurrence length would follow in 
the limit of continuous sampling. However, again we do not reach the regime and the 
theoretical limiting values are strongly dependent on the cut-off k, where the capillary 
wave theory breaks down. 



9 Dependence of p^ih) on n and h 

We plot in Fig. |9] the experimental curves for p^ih) for a large number of n for 
statepoint IV. A noteworthy point is that only for rather large values of n the decay 
with n (time) becomes exponential. The scatter in the data is modest, even for large 
n corresponding to large times t. Thus the experiment provides a host of detailed 
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Figure 10: The temporal Pn{h) as function of h for a number of n for statepoint IV. 
The drawn hues are the theoretical values. 



information on the statistics of the fluctuations in a wide time range. 

Another way of plotting the data is to select one value of n and plot Pnih) as 
function of h. Fig. 10 shows the experimental data for statepoint IV for Pn{h). This 
way of presenting the data facilitates the comparison with the theory. The drawn 
lines are the theoretical values as given by (21). We reiterate that the input in the 
calculations is a set of + 2 experimental values of g{0,t). The agreement between 
theory and experiment is remarkable for these detailed data. Statepoint II gives similar 
results with a slight asymmetry between up and down, of the same type as deviations 
in the spatial correlation functions, shown in Fig. |5} 

The data for p~{h) have been independently collected. The symmetry (25) is very 
well obeyed, such that there is no point in showing these data separately. 



10 Discussion 

The above given analysis of the statistics of interface fluctuations naturally falls into 
two parts, in which the height-height correlation function g{x,t) plays a pivotal role. 
The first part concerns the connection between g{x, t) and the state parameters such as 
Ap, rj and 7. The second part is the determination of the multiple correlation-functions 



G{hi,--- ,tn) from g{x,t) through (45). In the first part we have used the data for 
g{x,t) to determine the state parameters. Even though the derivation of the structure 
of g{x,t) in space and in time is quite different, the behavior is remarkably similar if 



the space and time variables are properly scaled (see (10) and (40)). 

The second part has been our main concern. We used the measured g{x,t) as in- 
put, providing all the necessary information on the statistics of the snapshots. The 
advantage of splitting the problem into these two parts is that the second part is not 
confounded by errors in the first. The only assumption in the theory is the use of the 
Langevin equation for the effect of the thermal (white) noise. The best justification 
for this procedure is a posteriori through its consequences. In view of the success- 
ful agreement with the experimental results, the assumption appears to be very well 
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fulfilled. 

Experimentally the capillary waves are disentangled from the structure of the in- 
trinsic interface. Most amazing is that detailed correlations in capillary waves can be 
determined with high accuracy. The statistics of the temporal dependence is generally 
better than that of the spatial behavior. We have chosen only a limited set of obtain- 
able correlation functions in order to compare them with the theoretical calculations. 
Experimentally it is easy to collect data for practically any interesting n. In Fig. |5]and 
Fig. [9] we show the distribution Pnih) as function of n for the values h = —1, and 1. 

There is a simplifying aspect in the fact that experimentally only sequences of 
finite time intervals can be measured. So one does not know what the interface does 
in between two snapshots. But this is precisely the reason why it suffices to calculate 
the correlation functions defined in Section |4j Here also one does not specify the 
evolution in between two snapshots. For instance, the key quantity q~^~{h,g) for the 
residence and waiting times, involves the crossing of the level h by the interface. But 
it does not say that it may cross it only once! Any odd number of crossings is possible. 
In Fig. [7| where we compare {h,g) with experiment, large time intervals feature 
(corresponding to small g) and in these large time intervals crossings are frequently 
taking place. Also the hills and valleys of length n for which the distribution is given 
in Fig. |5]and Fig. [T0| may be interrupted by opposite values in between snapshots. The 
charm of the comparison is that both theory and experiment allow these possibilities. 

In this paper we have restricted ourselves to sequences of height measurements 
at the same time or at the same position. The general result (45) shows that one 
could equally well correlate snapshots at combinations of times and positions and 
do a similar statistical analysis. The only point that matters is the height-height 
correlation gi j between the events i and j. Also one does not need to worry in how 
much the measurements refer to a single point or to an area of finite size. These more 
collective variables are also linear combinations of the basic variables /ik and therefore 
also gaussian randomly distributed. Then taking the measured correlations between 
the more general variables as input, leads to exactly the same analysis as given here 
for point variables. 
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A The height-height correlation integral 

In this appendix we discuss some properties of the height-height correlation function 
g{x,t). We start with the equal time function g{x,0). Using the scaled integration 
variable y = the integral for g{x, 0) leads to 

log(l + K^) Jo 1 + 
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where k = kmaxi is the cut-off. Sending this value to oo yields the modified Bessel- 
function Kq 

ydy:^^^Mf = Ko{x/0 (49) 
For finite k we can make a short distance expansion, reading 



ydy 



1 + 2/2 



log(l + ^^) + 



Matching the small argument expansion of Kq{x/^) 

Ko{x/0 = log(x/20 - 0.57721 + 



+ ■■■ (50) 



(51) 



with (50) gives for the point where they cross the approximate value 



X 



(52) 



The behavior in the time direction is remarkably similar to the spatial direction, 
although the integral for g{0,t) looks quite different. We write 



9{0,t) 

with H{z] k) as the integral 

H{z;k) = 



log(l + /s:2 



" ydy 
1 + r 



exp[-z{y + y'^)/2]. 



The first point is to prove relation (|40|), which we do by showing that 

-K^(z). 



dH{z, oo) 



dz 



(53) 



(54) 



(55) 



and checking the asymptotic expansion of (40) for large z. The advantage of (55) is 
that 

dH{z; oo) 



dz 



1 f°° 

- dyexp[-z{y + y-^)/2] 



(56) 



is a simpler integral than (54). Then take x = lay as integration variable, which turns 
(56) into 



dH{z] oo) 1 
dz ~ ~2 



dx exp(— 2; cosh x). 



(57) 



Splitting the integral into pieces from —00 to and from to 00 and changing in the 
first part from x to —x, yields the relation 



dH{z; 00) 
dz 



dx cosh X exp ( — z cosh x) . 



(5^ 



The integral is a representation of the function Ki{z) [16 



In order to show that no constant is lost in going from (40) to (55) one can check 
the asymptotic expansion of (40) for large z, which follows from an expansion around 
the slowest mode for ?/ = 1: 



y + y-' = 2 + {y-iy + 



(59) 
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and replacing the integral by a full gaussian around y = I. Then one gets 



H(z; oo) ^ — y diy- 1) 



exp[-(i/-l)V2] = e-^(^)'^', (60) 



which matches the asymptotic expansion of Kq{z). 

The final point is the expansion for small times t/tc for a finite value of k. We 
expand the exponential 

exp[-{y + y-')z/2] = 1 - (y + y-')z/2 + ■■■ (61) 



and insert this expansion into the integral (54). Then we find for H{z; k 



H{z-K) = hn{l + K^)-'^z + --- (62) 

Note that the next term in this expansion leads to a logarithmically divergent integral 
at the small y side. Thus the next term is not of the order z"^ but of the order z"^ In z. 
The finite k integral stays finite in contrast to Kq{z) which diverges for 2; = 0. With 



the expansion (51) and (62) we find for the point where the finite-^ curve starts to 
deviate from the Kq[z): 

t ~ tjK. (63) 

One obtains the rough estimate z ~ 1/k for this matching point, by looking to the 
value of the exponential at the upper boundary, which is exp[— z(k + /t~^)/2]. For 
larger values the boundary value starts to vanish and extending the integral to infinity 
leads to a small error. For smaller values of this z, the exponent of the exponential 



becomes smaller than 1 at the boundary and the integrand of (54) has not yet died out 



at y = K. Then deviations from the infinite domain start to show up. 
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